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Abstract 

We present two continuous symmetry reduction methods for reducing high- 
dimensional dissipative flows to local return maps. In the Hilbert polynomial 
basis approach, the equivariant dynamics is rewritten in terms of invariant coor- 
dinates. In the method of moving frames (or method of slices) the state space is 
sliced locally in such a way that each group orbit of symmetry-equivalent points 
is represented by a single point. In either approach, numerical computations 
can be performed in the original state-space representation, and the solutions 
are then projected onto the symmetry-reduced state space. The two methods 
are illustrated by reduction of the complex Lorenz system, a 5-dimensional dis- 
sipative flow with rotational symmetry. While the Hilbert polynomial basis 
approach appears unfeasible for high-dimensional flows, symmetry reduction by 
the method of moving frames offers hope. 
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1. Introduction 

In his seminal paper, E. Lorenz [1] reduced the continuous time and dis- 
crete spatial symmetries of the 3-dimensional Lorenz equations, resulting in a 
1-dimensional return map that yields deep insights [2] into the nature of chaos 
in this flow. For strongly contracting, low-dimensional flows, Gilmore, Lefranc 
and Letellier [3, 4] systematized construction of such discrete time return maps, 
through use of topological templates, Poincare sections (to reduce the contin- 
uous time invariance) and invariant polynomial bases (to reduce the spatial 
symmetries). They showed that in presence of spatial symmetries one has to 
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'quotient' the symmetry and replace the dynamics by a physicahy equivalent 
reduced, desymmetrized flow, in which each family of symmetry-related states 
is replaced by a single representative. This approach leads to symbolic dynam- 
ics and labeling of all periodic orbits up to a given topological period. Periodic 
orbit theory can then yield accurate estimates of long-time dynamical averages, 
such as Lyapunov exponents and escape rates [5]. 

In a series of papers Cvitanovic, Putkaradze, Christiansen and Lan [6-11] 
showed that effectively low-dimensional return maps can be constructed for 
high-dimensional (formally infinite dimensional) flows described by dissipative 
partial differential equations (PDEs) such as the Kuramoto-Sivashinsky equa- 
tion (KS). Such flows have state-space topology vastly more complicated than 
the Lorenz flow, and collections of local Poincare sections together with maps 
from a section to a section are required to capture all of the important asymp- 
totic dynamics. These KS studies were facilitated by a restriction to the flow-in- 
variant subspace of odd solutions, but at a price: elimination of the translational 
symmetry of the KS system and with it physically important phenomena, such 
as traveling waves. Traveling (or relative) unstable coherent solutions are ubiq- 
uitous and play a key role in organization of turbulent hydrodynamic flows, 
as pointed already in 1982 by Rand [12], and confirmed both by simulations 
and experimentation [13-18]. For KS [19, 20], and even for a relatively low- 
dimensional flow such as the complex Lorenz equations [21, 22] used as an ex- 
ample here, with the simplest possible continuous (rotational) spatial symmetry, 
the symmetry-induced drifts obscure the underlying hyperbolic dynamics. 

The question that we address here is how one can construct suitable return 
maps for arbitrarily high-dimensional but strongly dissipative flows in presence 
of continuous symmetries. Our exposition is based in part on refs. [5, 20, 23]. 
The reader is referred to Golubitsky and Stewart [24], Hoyle [25], Olver [26], 
Bredon [27], and Krupa [28] for more depth and rigor than would be wise to 
wade into here. 

In sect. 2 we review the basic notions of symmetry in dynamics. Sect. 2.1 
introduces the S0(2) equivariant complex Lorenz equations (CLE), a 5-dimen- 
sional set of ODEs that we use throughout the paper to illustrate the strengths 
and drawbacks of different symmetry reduction methods. In sect. 3 we describe 
important classes of solutions and their symmetries: equilibria, relative equilib- 
ria, periodic and relative periodic orbits, and use them to motivate the need for 
symmetry reduction. 

In sect. 4 we describe the problem of symmetry reduction. The action of 
a symmetry group endows the state space with the structure of a union of 
group orbits, each group orbit an equivalence class. The goal of symmetry 
reduction is replace each group orbit by a unique point a lower-dimensional 
reduced state space. In sect. 5 we briefly review the standard approach to spatial 
symmetry reduction, projection to a Hilbert basis, and explain why wc find it 
impracticable. In sect. 6 we review the method of moving frames, a direct 
and efficient method for computing symmetry-invariant bases that goes back to 
Cartan, and in sect. 6.1 we apply the method to the complex Lorenz equations. 
The method maps all solutions to a slice, a submanifold of state space that 
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plays a role for group orbits akin to the role Poincare sections play in reducing 
continuous time invariance. In contrast to the Hilbert basis approach, slices are 
local, with a generic trajectory within a slice bound to encounter singularities, 
and more than one slice might be needed to capture the flow globally. In sect. 6.2 
we show that a single local slice can suffice for the purpose of reducing the 
complex Lorenz equations flow to a return map. In sect. 7 we recast the method 
of moving frames into the equivalent, differential method of slices, with time 
integration restricted to a slice fixed by a given state space point. 

2. Symmetries of dynamical systems 

Here we arc interested in the role continuous symmetries play in dynamics. 
The methods we develop are in principle applicable to to translational and rota- 
tional symmetries of ODEs and PDEs, described by compact or noncompact Lie 
groups. We have in mind applications to PDEs such as Kuramoto-Sivashinsky 
and plane Couette flow which exhibit translational symmetries in either infinite 
or periodic domains. In the former case the group of symmetries is Euclidean 
and noncompact, and in the latter case orthogonal and compact. In numeri- 
cal computations the periodic setting is usually considered and, though Fourier 
analysis, a translation is represented by the action of the 1-parameter Lie S0(2) 
group on its linearly irreducible subspaces, the Fourier modes. Through trun- 
cation (for example spectral discretization), PDEs arc transformed to high- but 
finite-dimensional systems of ODEs. The key concepts will thus be illustrated 
by a specific ODE example, the S0(2) group acting on a five-dimensional state 
space, linearly decomposable into a direct sum of irreducible subspaces of SO (2). 

Consider a system of ODEs of the form 

X = v{x) (1) 

with V a smooth vector field and x ^ Ai cM.'^. 
A linear action g is a symmetry of (1) if 

v{gx)^gv{x) (2) 

for aU a; e M'^. One says that v commutes with g or that v is g-equivariant. 
When V commutes with the set of group elements g G G, the vector field v is 
said to be G-equivariant. The group G is said to be a symmetry of dynamics if 
for every solution x{t) = f^{x), g x{t) is also a solution. The finite time fiow 
/"^ (gx) through gx then satisfies the equivariance condition 

rigx)^gr{x). (3) 

In physics literature the term invariant is most commonly used; for example, in 
Hamiltonian systems a symmetry is manifested as invariance of the Hamiltonian 
under the symmetry group action. 

An clement of a compact Lie group continuously connected to identity can 
be written as 

g{d)^e''^, 0-T-;^0,T„ a- 1,2,... ,iV, (4) 
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where • T is a Lie algebra element, and 6a are the parameters of the trans- 
formation. Repeated indices are summed throughout this chapter, and the dot 
product refers to a sum over Lie algebra generators. The Euclidian product of 
two vectors x, y is indicated by x-transpose times y, i.e., x^y = Xiyi. Finite 
transformations exp(0 • T) are generated by sequences of infinitesimal steps of 
form 

g{Se)c^l + Sd -T, S9eR^, |(56i| < 1 , (5) 

where Ta, the generators of infinitesimal transformations, are a set of N linearly 
independent [dxd] anti-hermitian matrices, (T^)^ = — T^, acting linearly on 
the d-dimensional state space Ai. For G C 0{n) the generators can always be 
brought to real, antisymmetric form T"^ = — T. The fiow induced by the action 
of the group on the state space point x is given by the set of tangent fields 

taix)i = {Ta)ijXj . (6) 

These tangent fields are always normal to the 'radial' vector x, 

x^taix)^0. (7) 

For an infinitesimal transformation (5) the G-equivariance condition (2) be- 
comes 

dv 

v(x) ~ (l-e-T)v(x + 9-Tx) =v(x)~e-Tv(x) + —e-Tx. 

dx 

Thus the infinitesimal. Lie algebra G-equivariance condition is 

taiv)-A{x)taix)=0, (8) 

where A = dv/dx is the stability matrix. The left-hand side. 



( Ta-^{Tax) ] v{y) 



(9) 



is known as the Lie derivative of the dynamical flow field v along the direction 
of the infinitesimal group- rotation induced flow ta{x) = TqX. The equivariance 
condition (8) states that the two flows, one induced by the dynamical vector fleld 
V, and the other by the group tangent field t, commute if their Lie derivatives 
(or the Lie brackets or Poisson brackets) vanish. 

Any representation of a compact Lie group G is fully reducible, and invariant 
tensors constructed by contractions of are useful for identifying irreducible 
representations. The simplest such invariant is bilinear, 

• T = ^ G^"^ , (10) 

a 

(a) 

where Gj is the quadratic Casimir for irreducible representation labeled a, 
and l*-"-* is the identity on the a-irreducible subspace, elsewhere. The dot 
product of two tangent fields is thus a sum weighted by Casimirs, 

t(x)^.t(x') = EC^"^^^'^?4- (11) 
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If X is not invariant (fixed under group actions), t{x)'^ ■ t{x) is strictly positive. 
t{x)'^ ■ t{x'), however, can take either sign, or even vanish. 

2.1. An example: Complex Lorenz equations 

Consider a complex generalization of Lorenz equations, 

X ~ —ax + cry , V ~ {p ^ z)x — ay 

z = {xy* +x*y)/2~bz, (12) 

where x, y are complex variables, z is real, while the parameters cr, h are real and 
p = pi + ip2, a = 1 — ie are complex. Recast in real variables, x = xi + ix2 , y = 
yi + iy2 this is a set of five coupled ODEs 

XI = -axi + ayi , ±2 = -0X2 + cr?/2 
2/1 = {p\ - z)xi - P2X2 - 2/1 - 62/2 

2/2 = P2X1 + {pi - z)x2 + eyi - y2 

z = -bz + xiyi + X2y2 ■ (13) 

In all numerical examples that follow, the parameters will be set to pi ~ 28, p2 = 
0, b ~ 8/3, cr ~ 10, e = 1/10, unless explicitly stated otherwise. Why worry 





Figure 1: State space portrait of complex Lorenz flow. Plotted are a generic chaotic trajectory 
(blue), the Eq equilibrium, a representative of its unstable manifold (green), the Qi relative 
equilibrium (red), its unstable manifold (brown), and one repeat of the 01 relative periodic 
orbit (black) (color online). 

about continuous symmetries? The visualization in figure 1 of typical long-time 
dynamics of complex Lorenz flow suffices to illustrate the effect a continuous 
symmetry has on dynamics. A generic trajectory slowly drifts along the direc- 
tion of continuous symmetry while tracing a Lorenz-butterfly like attractor. It 
is a mess. 

The complex Lorenz equations are a dynamical system with a continuous 
(but no discrete) symmetry, equivariant under the one-parameter rotation group 
U(l) ^ S0(2) acting by 

(x, y, z) ^ {e^'x, z) , e [0, 2n] . (14) 
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Alternatively, substituting the Lie algebra generator 



/ -1 \ 

1 

-1 

10 

V y 

acting on a 5-dimensional space (13) into (4) yields the 
finite angle action (14) of S0(2) 



(15) 
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(16) 



We see that the linear action of SO (2) on the state space of the complex Lorenz 
equations decomposes into the m = G-invariant subspace (z-axis) and the 
m ~ 1 subspace of multiplicity 2. 

The generator T is anti-symmetric, T-^ = — T, and the group is compact, 
its elements parametrized by 9 mod 2tt. Locally, at a; £ M, the infinitesi- 
mal action of the group is given by the group tangent field t{x) = Tx = 
{—X2,xi,—y2,yi,0). In other words, the flow induced by the group action is 
normal to the radial direction in the {xi, X2) and (yi, 2/2) planes, while the z-axis 
is left invariant. 

The equivariance of the complex Lorenz flow under S0(2) rotations (16) 
can be verified by substituting the Lie algebra generator (15) and the stability 
matrix for complex Lorenz flow (13), 





Pi - z 

P2 

yi 





— (T 

Pi - Z 

y2 



a 


-1 

e 





a 






— e —x\ 

-1 -~X2 

X2 -b J 



(17) 



into the equivariance condition (8). For the parameter values (13) the flow is 
strongly volume contracting. 



tiA = ^ Ai(a;,t) = -6 - 2(cr + 1) = -24-2/3. 



(18) 



The complex Lorenz equations (12) were introduced by Gibbon and McGuin- 
ness [21, 22] as a low-dimensional model of baroclinic instability in the atmo- 
sphere. Zeghlache and Mandel [29] and Ning and Haken [30] have shown that 
equations isomorphic to the complex Lorenz equations, with e + p2 ~ 0, also 
appear as a truncation of Maxwell-Bloch equations describing a single mode. 
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detuned, ring laser. The choice e + /92 = is degenerate (see (22)) in the sense 
that it leads to non-generic bifurcations. We follow Bakasov and Abraham [31] 
who set p2 =0 and e 7^ to describe detuned ring lasers. 

Here, however, we are not interested in the physical applications of these 
equations; rather, we study them as a simple example of a dynamical system 
with continuous (but no discrete) symmetries, with a view of testing methods of 
reducing the dynamics to a lower-dimensional reduced state space. We investi- 
gate various ways of quotienting its S0(2) symmetry, and reducing the dynamics 
to a 4-dimensional reduced state space. As we shall show, the dynamics has a 
nice 'stretch & fold' action, but that is totally masked by the continuous sym- 
metry drifts. We shall not rest until we attain the simplicity of figure 6, and 
the bliss of 1-dimensional return map of figure 4. 

3. Symmetries of solutions 

In order to explore the implications of equivariance on solutions of dynamical 
equations, we start by examining the way a compact Lie group acts on a state 
space M. The group orbit or G-orhit of the point x ^ M. \s the set 

M^^{gx\geG} (19) 

of all state space points into which x is mapped under the action of G. The 
symmetry Gx {isotropy or stabilizer group) of a state space point x is the largest 
subgroup of G 

Gx ^ {g e G : gx ^ x} (20) 

that leaves x fixed. The symmetry Gx of a set Mx € is the largest subgroup 
of G that leaves M.x invariant as a set: 

Gx = {g ■■ gMx = Mx} ■ 

If Gp is a symmetry, intrinsic properties of a solution Mp (such as equilibrium 
or a cycle stability eigenvalues, period, Floquet multipliers) evaluated anywhere 
along its Gp-orbit are the same. A symmetry thus reduces the number of in- 
equivalent solutions. So we also need to describe the symmetry of a solution, as 
opposed to (3), the symmetry of the system. 

The fixed-point subspace Fix (H) of a subgroup H d G is the subspace of Ai 
containing all fixed points of H: 

Fix (H) ^ {x e M, g e H \ gx = x} . 

The physical importance of fixed-point subspaces lies in the fact that they are 
invariant under G-equivariant dynamics [24], 

r (Fix (H)) C Fix (H) 

and thus flow invariant for all times t. Therefore if x{t) is a solution of an 
equivariant ODE, then its symmetry Gx{t) = Gx(q) is preserved for all times. 
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Figure 2: (a) A relative equilibrium orbit starts at some point x(0), with the dynamical flow 
field v{x) = c-t(x) pointing along the group tangent space. For the SO(2) symmetry depicted 
here, the flow traces out the group orbit of x{0) in time T= 2it/c. An equilibrium lives either 
in the Fix(G) subspace (0:3 axis in this sketch), or on a group orbit as the one depicted here, 
but with zero angular velocity c. In the latter case the circle (in general, TV-torus) depicts a 
continuous family of flxed equilibria, related only by the group action, (b) A relative periodic 
orbit starts out at x{0) with the dynamical v and group tangent t flows pointing in different 
directions, and returns to the group orbit of x{0) after time Tp at x{Tp) = gpx{0), a rotation 
of the initial point by gp. 

In contrast to equilibrium solutions that satisfy f^ix) ~ x, relative equilibria 
(or traveling waves) satisfy /"^(x) = g{T) x for any r, where the group has been 
reparameterized by time, 9 = 0{t). In a co-moving frame moving along the 
group orbit with velocity v{x) = c • t{x), the relative equilibrium appears as an 
equilibrium. Here t is the group tangent field (6). 

A relative periodic orbit is an orbit Aip for which the initial point exactly 
recurs 

Xp{Q) = gpXp{Tp) , Xp{T)eMp, (21) 

at a fixed relative period Tp, but shifted by a fixed group action gp which brings 
the endpoint Xp(Tp) back into the initial point Xp{0), see figure 2(b). The 
group action gp = gp{9) parameters 9p = {di,02, ■ ■ -d^) wiU be referred to as 
phases, or shifts. For dynamical systems with only continuous (no discrete) 
symmetries, the parameters {t,9i,--- ,0^} are real numbers, the ratios tt/Bj 
are almost never rational, and the likelihood of closing into a periodic orbit is 
zero. Thus the trajectory of a relative periodic orbit generically sweeps out the 
group orbit ergodically. 

A relative periodic orbit is periodic in its mean velocity Cp = 9p/Tp co- 
rotating frame, figure 3, but in the stationary frame its trajectory is quasiperi- 
odic. A co-moving frame is helpful in visualizing a single 'relative' orbit, but 
useless for viewing collections of orbits, as each one drifts with its own group 
velocity. A simultaneous visualization of all relative periodic orbits as periodic 
orbits can be attained only by symmetry reduction, to be undertaken in sects. 5 
and 6. 

Relative equilibria and relative periodic orbits are the hallmark of systems 
with continuous symmetry. Amusingly, in this extension of periodic orbit theory 
from unstable 1-dimensional closed periodic orbits to unstable (A^-l-l)-dimens- 
ional compact manifolds Aip invariant under continuous symmetries, there are 
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Figure 3: A relative periodic orbit of the Kuramoto-Sivashinsky flow, traced for four periods 
Tp and projected on (a) a stationary state space coordinate frame {di, D2, V3}; (b) a co-moving 
{^'l , , 'I'a } coordinate frame, moving with the mean velocity Cp = 6p/Tp. (Prom ref. [19].) 

either no or proportionally few periodic orbits. In presence of a continuous and 
no discrete symmetry, likelihood of finding a periodic orbit is zero. Relative 
periodic orbits are almost never eventually periodic, i.e., they almost never lie 
on periodic trajectories in the full state space, so looking for periodic orbits in 
systems with only continuous symmetries is a fool's errand. 

A historical note. Relative equilibria and relative periodic orbits are related 
to equilibria and periodic orbits of dynamics reduced by the symmetries. They 
appear in many physical situations, such as motion of rigid bodies, gravita- 
tional A^-body problems, molecules, nonlinear waves, spiralling patterns and 
turbulence. According to Cushman, Bates [32] and Yoder [33], C. Huygens [34] 
understood the relative equilibria of a spherical pendulum many years before 
publishing them in 1673. A reduction of the translation symmetry was obtained 
by Jacobi (for a modern, symplectic implementation, see Laskar et al. [35]). Ac- 
cording to Chenciner [36], the first attempt to find (relative) periodic solutions 
of the TV-body problem was the 1896 short note by Poincare [37], in the context 
of the 3-body problem. Relative equilibria of the TV-body problem (known in 
this context as Lagrange points, stationary in the co-rotating frame) are cir- 
cular motions in the inertial frame, and relative periodic orbits correspond to 
quasiperiodic motions in the inertial frame. Relative equilibria that exist in a 
rotating frame are called central configurations. For relative periodic orbits in 
celestial mechanics see also ref. [38]. A striking application of relative periodic 
orbits has been the discovery of 'choreographies' of TV-body problems [39-41]. 

The modern story on equivariance and dynamical systems starts perhaps 
with M. Field [42], and on bifurcations in presence of symmetries with Ru- 
elle [43]. Ruelle proves that the stability matrix/ Jacobian matrix evaluated at 
an equilibrium/fixed point x € M.g decomposes into linear irreducible represen- 
tations of G, and that stable/unstable manifold continuations of its eigenvectors 
inherit their symmetry properties, and shows that an equilibrium can bifurcate 
to a rotationally invariant periodic orbit (i.e., relative equilibrium). 
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3.1. An example: Solutions of the complex Lorenz equations 

In the case of the complex Lorenz equations the origin Eq is an equihbrium 
of (12) for any value of the parameters. It is stable for < pi < pic and unstable 
for pic < pi, where [22] 



Pic 



l + {e + P2){e^ap2)/{cr + l)^ 



At the bifurcation [43] a pair of eigenvalues crosses the imaginary axis with 
imaginary part 

Uc = (7{e + p2)/i<J + l), (22) 

and a relative equilibrium Qi with constant angular velocity Wc is born. For 
ujc = the relative equilibrium degenerates to an SO(2)-orbit of equilibria. As 
the existence of a relative equilibrium in a system with S0(2) symmetry is the 
generic situation, we follow ref. [31] and set p2 = and e 7^ 0. 

To find the location of the relative equilibrium it is convenient to work in 
polar coordinates 

{xi,X2, yi,y2,z) ^ (ri cos 6*1, ri sin^i, r2 cos 6*2, r2 sin 02, 2) , (23) 

where > , r2 > 0. The complex Lorenz equations (12) take the form 

/ -cr (ri - r2 cos6') \ 



Oi 
r2 
O2 
V i J 



V 



ri 

-r2 + ri {{pi — z) cos 9 — p2 sin ( 
e + 77 ((Pi - z) sin6' + p2 cos 6') 

—bz + rir2 cos 9 



J 



For rotationally invariant flows the dynamics depends only on the relative angle 
9 = 9i — 92 (which is why one speaks of 'relative' equilibria). This observation 
enables us to recast the complex Lorenz equations in the 4-dimensional reduced 
state space: 



[ ri\ 



\ i J 



( 



\ 



—a (ri — r2 cos 6') 
-"^2 + {pi — 'Ari cost 



ar^ + ipi-z)^^jsm( 
—bz + rir2 cos 9 



(24) 



where we have set p2 = 0- The full 5-dimensional evolution can be regained by 
integrating the driven reconstruction equation for the mean angular velocity: 



?i + 6*2 = e - [ar2/ri - {pi - z)ri/r2] sint 



(25) 



In general 9i and ^2 change in time, but for the relative equilibria the difference 
between them is constant. The condition for a relative equilibrium is that all 
time derivatives in (24) vanish, while ^1 = 02 7^ (if 0i = 02 = we have a 
group orbit of equilibria instead). The relative equilibrium Qi is given by 

(ri,r2,0,z) = ( {Pi - d), ^/bd {pi - d), cos-\l/Vd),pi - d) , (26) 
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where d = 1 + e^/ {a + 1)^, and its angular velocity is 

e, = ae/{a + 1) , (27) 

with period Tq^ — 27r(a + l)/ae. For the parameter values (13), the relative 
equilibrium is at 

XQ, = {ri,r2,e,z) = (8.48527,8.48562,0.00909,26.9999), (28) 

rotating with the period Tqj = 69.1150. 

As pi is increased, a secondary bifurcation from Qi results in a relative pe- 
riodic orbit (21), or, more precisely, in the quasiperiodic 2-frequency modulated 
traveling wave [28]. With further increase in pi the dynamics turns chaotic, 
with an infinity of unstable relative periodic orbits. Once symmetry reduced 
maps are constructed (see figure 4(b) below), a large number of these can be 
computed by methods described elsewhere [19, 20]. Calculation of the Qi sta- 
bility eigenvalues for the parameter values (13) (sec rcf. [20] for a calculation of 
stability of relative equilibria in cquivariant variables) yields a weakly unstable 
spiral-out equilibrium 

(Ai,2, A3, A4) = (0.0938179 ± 10.1945i, -11.0009,-13.8534) . (29) 

The role of the above exact invariant solutions is illustrated by the portrait 
of complex Lorenz fiow state space in figure 1, with the relative equilibrium 
Qi and three repetitions of the 01 relative periodic orbit superimposed over 
a generic chaotic orbit. Repeats of 01 trace out a torus crgodically, so in a 
system with a 1-dimensional continuous symmetry the organizational blocks of 
a strange attractor are circles (relative equilibria) instead of points (equilibria), 
and partially hyperbolic tori (relative periodic orbits) instead of closed loops 
(periodic orbits) . It is difficult to understand the geometry of the flow by looking 
at such tori. 

The large imaginary part of Ai in (29) implies that the simulation has to be 
run up to time of order of at least 70 for the strange attractor in figure 1 to start 
filling in. Dynamics is organized by the interplay of the stable and unstable 
manifolds of equilibrium Eq and relative equilibrium Qi, but the symmetry- 
induced drift along the direction of rotation blurs the picture and the notion of 
recurrence becomes relative. In what follows, it is this confusing situation (as 
well as the theoretical fact [44] that dynamical zcta functions have their support 
on relative periodic orbits) that motivates the search for effective methods to 
project the dynamics onto a reduced state space. 

4. Symmetry reduction 

The action of a symmetry group G on endows the state space with the 
structure of a union of group orbits, each group orbit an equivalence class. 
The goal of symmetry reduction is the identification of a unique point as the 
representative of a group orbit, and the replacement of the original state space 
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by the space of such points, the reduced state space. In the hterature this space 
is alternatively called desymmetrized state space, symmetry-reduced space, orbit 
space, or quotient space Ai/G because symmetry has been 'divided out.' The 
symmetry group G of equivariant dynamics acts trivially in the reduced state 
space, and the resulting dynamical system, called by Gilmore and Lettelier [4] 
the image, is symmetry invariant, in the sense that its symmetry group is the 
identity. Reduced state space is in general not a manifold but rather a union of 
manifolds of different dimensions [45] . 

In sect. 5 we briefly review one of the standard tools by which spatial sym- 
metry reduction can be achieved: projection to a Hilbert basis, and show in 
sect. 5.1 how it works for complex Lorenz flow. A wonderful symmetry reduc- 
tion tool for low-dimensional flows, the Hilbert basis approach turns out to be 
too cumbersome to be applicable to high-dimensional flows. Next we describe 
the method of moving frames (sect. 6) and apply it to the complex Lorenz flow 
example to illustrate the form of a general linear slice (sect. 6.1), show how 
the method enables us to explicitly compute G-invariant coordinates, and re- 
late these to the Hilbert invariant polynomial basis (sect. 6.2). Then we discus 
difii'erent choices of slice- fixing points (sect. 6.3), and the associated singular 
sets. Since rotations commute with time integration, one can start with a point 
on the slice, integrate for short time and then rotate the end point back into 
the slice. In sect. 7 the limit of infinitesimal time steps yields the equivalent 
but differential formulation, the method of slices for which the flow is restricted 
to the reduced state space. In practice we find it more convenient to use the 
numerical code as given, and post-process the data by the method of moving 
frames, rather than rewriting the equations in the method of slices form. 

5. Hilbert polynomial bases 

In atomic physics and other low-dimensional physical problems with spatial 
symmetries, symmetry reduction is customarily implemented just as we did it in 
(23), by going to the natural coordinate system (polar, cylindrical, etc.). That 
works well for linear systems, but not so well for nonlinear flows, and some 
take pride in using no polar coordinates in symmetry reduction of Hamiltonian 
flows [32, 46]; note, for example, that these coordinate transformations introduce 
singularities in (24) at ri = r2 = 0. 

What are we really doing when redefining dynamics in terms of such invariant 
coordinates? We are recasting equivariant dynamics of {xi,X2, - • •) coordinates 
in terms of rotationally invariant lengths (ri = {xf + x^Y , ■ ■ ■), volumes and 
other invariant quantities. Physical laws have the same form in all coordinate 
frames, so they arc often formulated in terms of functions (Hamiltonians, La- 
grangians, • • • ) that are invariant under a given set of symmetries. Given a 
symmetry, what is the most general functional form of such law? The general 
problem of symmetry reduction in this sense was elegantly solved nearly a cen- 
tury ago. According to the Hilbert- Weyl theorem, for a compact group G there 
exists a finite G-invariant Hilbert polynomial basis M2, . . . , itm}i m > d, 
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such that any G-invariant polynomial can be written as a multinomial 



h{x) = p{ui{x),U2{x), . . .,u„,{x)) 



X e M. 



(30) 



The Gilmore and Lettelier monograph [4] offers a clear, detailed and user friendly 
discussion of symmetry reduction by means of invariant polynomial bases (do 
not look for Hilbert in the index, though). The invariant dynamical equations 
follow from the equivariant ones by chain rule 



upon substitution {xi, X2, • • • , Xd} — > {ui, W2, • • • , u„i}. One can either rewrite 
the dynamics in this basis, or one can simply plot the 'image' of solutions com- 
puted in the original, equivariant basis in terms of these invariant polynomials. 

Unfortunately, while the idea is elegant, an explicit construction of G-invari- 
ant basis can in practice be a daunting undertaking. The set of m > d invariant 
polynomials {ui, M2, ■ ■ ■ , Um} is not unique, and while these polynomials are lin- 
early independent, they are functionally dependent through m — d+N nonlinear 
relations called syzygies. Their determination becomes quickly computationally 
prohibitive as the dimension of the system and/or group increases [45, 47], and 
in practice computations are confined to dimensions less than ten. As our goal 
is to quotient continuous symmetries of high-dimensional flows, as high as 10^- 
10^ coupled ODEs arising from truncations of the Kuramoto-Sivashinsky and 
Navier-Stokes flows, reduction by the method of Hilbert basis is at present not 
a feasible option. 

Nevertheless, as symmetry reduction of moderate-dimension flows by the 
method of invariant polynomials offers a clean benchmark for other approaches 
to symmetry reduction, we start by showing how it works for complex Lorenz 
flow. 

5.1. An example: Complex Lorenz equations recast in Hilbert basis 

As the Hilbert basis approach turns out to be too cumbersome for our main 
goal, symmetry reduction of high-dimensional flows, we forgo here a systematic 
discussion of how to construct invariant polynomials bases. The pedagogical 
literature mostly focuses on discrete symmetry groups [4, 45, 47], while general 
algorithms are a domain of advanced algebraic geometry monographs. For the 
purpose at hand it suffices to use the Gilmore and Letellier [4, 48] invariant 
polynomial basis for the action (16). They apply it to the symmetry reduc- 
tion of the Zeghlache-Mandel system [29] , a flow much like the complex Lorenz 
equations. As it can be easily verified, the Hilbert basis 



(31) 



dXj 




2 I 2 

Ui = xiyi + X2y2 



U3 = xiy2 ~ X2yi , 



(32) 



U5 ^ z 
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is invariant under (16), the S0(2) action on a 5-diniensional state space. That 
imphes, in particular, that the image of the full state space relative equilibrium 
Qi group orbit of figure 1 is the equilibrium point in figure 4 (a), while the image 
of a relative periodic orbit, such as 01, is a periodic orbit. The five polynomials 
are linearly independent, but related through one syzygy, 

U1U2 — U3 — U4 = , (33) 

yielding a 4-dimensional 7M/SO(2) reduced state space, a symmetry-invariant 
representation of the 5-dimensional S0(2) equivariant dynamics. Vladimirov, 




Figure 4: (a) Invariant image of complex Lorenz flow, figure 1, projected onto the invariant 
polynomials basis (32). Note the unstable manifold connection from the equilibrium Eq at 
the origin to the strange attractor controlled by the rotation around the reduced state space 
image of relative equilibrium Qi; as for the Lorenz flow [1], the natural measure close to Eq 
is vanishingly small but non-zero, (b) The return map for the Poincare surface of section 
Ul = U4 for the complex Lorenz equations projected on invariant polynomials (32). The 
return map coordinate is the Euclidean length (36) along the Poincare section of the unstable 
manifold of Qi. 

Toronov and Dcrbov [49] use a different invariant polynomial basis to study 
bounding manifolds of the symmetry reduced complex Lorenz flow and its ho- 
moclinic bifurcations. 

Application of the chain rule (31) brings the equivariant complex Lorenz 
equations (13) to the invariant form (32): 



ill = 2 (T (U4 — Ul) , 

U2 = -2{U2 - P2U3 - {pi - U5)U4) , 

U3 = -(cr + 1) U3 + P2 "1 + eM4 , (34) 

U4 = -((T + 1) W4 + (Pl - M5) "1 + — e U3 , 

U5 = li4 — 5 li5 . 



As far as visualization goes, we need neither construct the invariant equations 
(34) nor integrate them. It suffices to integrate the original, unreduced flow of 
Figure 1, but plot the solution in the image space, i.e., the invariant, Hilbert 
polynomial coordinates Ui, as in figure 4(a). A minor drawback of the Hilbert 
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polynomial basis projections is that the folding mechanism is harder to view 
since the dynamics is squeezed near the z-axis. 



5.2. Symmetry reduced return map 

Successive trajectory intersections with a Poincare section^ a (d — l)-dim- 
ensional manifold or a set of manifolds V embedded in the d-dimensional state 
space A4, define the Poincare return map P{x), a (d — l)-dimensional map of 
form 

x' = P(a;) = r(^)(x), x',xeV. (35) 

Here the first return function t{x), sometimes referred to as the ceiling function., 
is the time of fiight to the next Poincare section for a trajectory starting at x. 
A good choice of the Poincare section manifold V in high-dimensional (d > 3) 
flows is basically a dark art. We chose a Poincare section which contains the 
z-axis and the relative equilibrium, here defined by the condition ui ~ 1*4. Even 
though in the complex Lorenz equations case we could use one of the Ui s as a 
coordinate to construct a return map, for high dimensional fiows a dynamically 
intrinsic parametrization is the only option. Following ref . [6] , we construct the 
first return map of figure 4 (b) using as a coordinate the Euclidean length along 
the intersection of the unstable manifold of Qi within the Poincare surface of 
section, measured from Qi. 

We begin by sprinkling evenly spaced points {a;^^^ , x*-^-* , • • • , a;*^^"^-* } between 
the relative equilibrium point xq^ = a;^°^ and the point x = x^'^\ along the Id 
intersection of Poincare section and unstable manifold continuation x''^'^ G W^^"-) 

of the unstable e^^^ eigenplane (we shall omit the eigendirection label (1) in what 
follows). Then the arclength in Euclidean metric from the relative equilibrium 
point XQ^ = x^°^ to X = x^^^ is given by 



lim 

N^oc 



\ 



N 



^(dxW) . (36) 



By definition /'^('^^(x) € M^^" , so /* (x) induces a Id map s(so, ''') = s(/'"(^°) (xq)). 

Turning points are points on the unstable manifold for which the local un- 
stable manifold curvature diverges for forward iterates of the map, i.e., points at 
which the manifold folds back onto itself arbitrarily sharply. The Id curve W,") 
starts out linear at xq^ , then gently curves until it folds back sharply at the 
'turning point' along a possible heteroclinic connection to Eg, and then nearly 
retraces itself. 

The trick is to figure out a good base segment to the nearest turning point 
L = [0,Sb\, and after the foldback assign to s{x,t) > Sf, the nearest point s on 
the base segment. Since here the stable manifold contraction is strong, the 2nd 
coordinate connecting .s{x,t) — > s can be neglected. 

Armed with this intrinsic curvilinear coordinate parametrization, we are now 
in a position to construct a 1-dimensional model of the dynamics on the non- 
wandering set. If x„ is the nth Poincare section of a trajectory in neighborhood 
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of Xq, and s„ is the corresponding curvilinear coordinate, then s„+i = /^"(sn) 
models the full state space dynamics Xn — >■ Xn+i- We approximate /(s„) by a 
smooth, continuous 1-dimensional map / : ig — > by taking x„ G Lq, and 
assigning to i„+i the nearest base segment point s„+i = s(a;„4-i). Thanks to 
the extreme contraction rate (18), the return map turns out to be unimodal for 
all practical purposes, so binary symbolic dynamics are easily constructed and 
admissible periodic orbits of the map up to desired length can be systemati- 
cally obtained. A multiple shooting routine [5] can then be used to determine 
the corresponding relative periodic orbits of the complex Lorenz equations to 
machine precision. 

6. The method of moving frames 

The method of moving frames, introduced by G. Darboux and systematized 
by E. Cartan [50], is interpreted by Fels and Olver [51, 52] as a map from a d- 
dimensional manifold to a A^-dimensional Lie group acting on it. Moving frames 
are then used to compute d—N functionally independent fundamental invariants 
for general group actions in relation to general equivalence problems. 'Funda- 
mental' here means that they can be used to generate all other invariants, and, 
in particular, they serve to distinguish group orbits in an open neighborhood of 
the slice point, i.e., two points lie on the same group orbit if and only if all fun- 
damental invariants agree. For an introduction to the method wc recommend 
Olver's pedagogical monograph [26]. Here wc emphasize the application of the 
method to dynamical symmetry reduction, and focus in particular on groups 
acting on spaces with the structure of a direct sum of irreducible subspaces, an 
application which is not to our knowledge explored in the literature. We are not 
concerned with the explicit determination of the fundamental invariants as in 
ref. [51, 52], except for illustrative purposes; instead we focus on implementation 
of a moving frame transformation as a numerically fast and efficient linear map 
of the full state space dynamics onto its desymmetrizcd, reduced state space 
projection. 

The main idea behind method of moving frames is that we can, at least 
locally, map each point along any solution x{t) to a unique representative x{t) 
of the associated group orbit equivalence class, by a suitable rotation 

x{t) ^ g{T)x{T) . (37) 

Equivariance implies the two points are equivalent. In the method of moving 
frames the reduced state space representative x of a group orbit equivalence 
class is picked by slicing across the group orbits by a fixed manifold. 

In the following it will be useful to introduce the notion of a slice, an (d—N)- 
dimensional submanifold Ai C A4 such that A4 intersects all group orbits in an 
open neighborhood otx' € A4 transversally and at most once. In other words, 
slice is the analogue of a Poincare section, but for group orbits. As is the case 
for the dynamical Poincare sections, in general a single slice does not suffice to 
intersect all group orbits of points in A4 . One can construct a local slice passing 
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X2-X2 




Figure 5: (a) Slice Ai is a hypcrplanc (38) passing through the slice- fixing point x' , and 
normal to the group tangent t' at x' . It intersects all group orbits (indicated by dotted lines 
here) in an open neighborhood of x' . The full state space trajectory x{t) and the reduced 
state space trajectory x{t) belong to the same group orbit and are equivalent up to a 

group rotation (?(t). (b) The method of moving frames for a flow SO(2)-cquivariant under 
(16) writh slice through x' = (0,1,0,0,0), group tangent t' = (—1,0,0,0,0). The orientation 
condition restricts the slice to half-hyperplane xi = 0, X2 > 0. A trajectory started on the 
slice at 3:(0) evolves to a state space point with a non-zero xi(ti). Compute angle 61 through 
slice condition (39). Rotate x(ri) counter-clockwise by 61 to x(ti) = g{9i) x{ti), so that the 
equivalent point on the circle lies on the slice, xi{ti) = 0. Repeat for all sample points x{Ti) 
along the trajectory. 

through any point x E Ai ii the group orbits of G have the same dimension, 
i.e. they are away from fixed-point subspaccs of continuous subgroups of G, see 
ref. [52] for details. 

The simplest slice condition defines a linear slice as a (d — Af)-dimension- 
al hyperplane Ai normal to the N group rotation tangents t'^ at point x' , see 
figure 5: 

{x~x'ft'^=0, t'^=ta{x') = Tax' , a = l,2,---,N. (38) 

As x''^t'i^ = by the antisymmetry of Tq, the slice condition (38) fixes 6 for a 
given x by orthogonality, 

= x'^t', = x'^g{0ft', , (39) 

where g'^ denotes the transpose of g. A group orbit will in general intersect a 
slice more than once (for example in the case of S0(2), two 7r-separated points), 
so we need to impose further conditions on the slice, in the form of either 
inequalities or orientation conditions, so as to ensure unique intersection. These 
restrictions are rather arbitrary, the only requirement being that M remains 
a connected manifold. We illustrate this point for the complex Lorenz flow 
example in sect. 6.1. 

For group orbits intersected by a slice, we can identify the unique group 
element g = g{x) that rotates x into the slice, gx = x € M. The map that 
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associates to a state space point x a Lie group action g[x) is called a moving 
frame. The method of moving frames can be thought of as a change of variables 
X = g^^{x)x , to a frame of reference for which the slice-fixing condition (39) is 
identically satisfied-hence the name 'moving frame.' 

The method of moving frames is a post-processing method; trajectories are 
computed in the full state space, then rotated into the slice whenever desired, 
with the slice condition easily implemented. The slice group tangent t' is a given 
vector, and rotation parameters 9 are determined numerically, by a Newton 
method, through the slice condition (39). For given 9, g(9) x is another vector, 
linear in x. 

A slice can be identified with M./G in an open neighborhood of x' . As is 
the case for the dynamical Poincare sections, a single slice does not suffice to 
reduce M ^ M/G globally as one cannot expect the group orbit of any point 
in M to intersect a given slice. 

How does one pick a slice point x'1 A generic point x' not in an invariant 
subspace (on the z axis of the complex Lorenz equations, for example) should 
suffice to fix a slice. The rules of thumb are much like the ones for picking 
Poincare sections. The intuitive idea is perhaps best visualized in the context 
of fluid flows. Suppose the flow exhibits an unstable coherent structure that is 
frequently visited at different spatial dispositions. One can fit a 'template' to 
one recurrence of such structure, and describe other recurrences as its transla- 
tions. A well-chosen slice point belongs to equivalence class that is dynamically 
important in this way (i.e., a group orbit). We discuss, in the context of our 
complex Lorenz equations example, several slice fixing choices in sects. 6.2 and 
6.3. 

A historical note. For the definition of slice see, for example, Chossat and 
Lauterbach [45] . Slices tend to be discussed in contexts much more difficult than 
our application - symplcctic groups, sections in absence of global charts, non- 
compact Lie groups. We follow ref. [53] in referring to a local group-orbit section 
as a slice. The usage goes back at least to Palais [54] in 1961 and Mastow [55] 
in 1957. Some [27, 52] refer to global group-orbit sections as cross-sections, a 
term that we would rather avoid, as it has an honest, well-established meaning 
in physics. Guillemin and Sternberg [56] define the cross-section, but emphasize 
that finding it is very rare: "existence of a global section is a very stringent 
condition on a group action. The notion of slice is weaker but has a much 
broader range of existence." 

6.1. An example: Moving frame for complex Lorenz equations 

In case of the complex Lorenz equations we can, due to equivariance, rotate 
any slice fixing point in (39) so that we have x'^ ~ 0. As only the group tangent 
direction matters, a slice that goes through point (0,^^21^15 ^2''^') equivalent 
to (0, l,y']^/x2, ^2/^21 0)1 so can specify the most general slice fixing point 
for the complex Lorenz equations by two numbers, 

5;' = (0, l,f;,f^, 0). (40) 
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The group orbit tangent then becomes t' = (—1, 0, —y'2, y'l, 0) and shce condi- 
tion (39) leads to 

= tan — — — . (41) 

X2 + VlVl + V2V2 

To ensure a unique intersection with the slice, we have to further restrict M. 
by choosing a representative out of the two group orbit points that intersect 
the slice. One can impose an orientation condition, for example choosing the 
point that is at minimum distance from x', or one can define the inverse tangent 
function tan~^(6/a) so that it distinguishes quadrants in the (a, h) plane. Either 
condition works equally well. 

We observe that (41) is undefined when 

+ 2/2^1 - y'lVi = , (42a) 
X2 + y'lVi + 2722/2 = , (42b) 

arc both satisfied. We will refer to this 3-dimensional linear subspace as the 
singular set of the moving frame associated with (41). Condition (42a) implies 
that point x is already on the slice, x'^x' = 0. Condition (42b) implies that 
the group tangent at point x is perpendicular to group tangent at slice fixing 
point, t[xyt' = —x'^x' = 0. The problem lies in the fact that the limit of (42) 
as we approach the singular set does not exist. For instance, consider points 
for which (42b) holds; as we approach the singularity with positive values of 
the numerator in (41) we have 9 = 7r/2, while for negative values 9 = — 7r/2. 
A TT-jump occurs as we cross the singular set. In general such crossings are 
expected to occur, since the singular set is not flow invariant. Even worse, the 
singularity distorts the way trajectories are mapped onto the slice, even if they 
merely approach it rather than cross it, as we will see in the next two examples. 

6.2. Irreducible subspace slice and explicit invariants 

We now show how a particular choice of the slice point enables us to express 
the transformation to invariant variables in a simple analytic form. Place the 
slice point in one of the linearly irreducible subspaccs of S0(2) action (16), 
for instance x' = (0,-1,0,0,0). The group-tangent at the slice point is then 
t' = (1, 0, 0, 0, 0) and the slice fixing condition is 

xi = xi cos 9 ~ X2 sin 9 ~Q . (43) 

The orientation condition restricts the slice to half-hyperplane xi = X2> 0. 
Solving (43) for the polar angle 9 in {xi,X2) we get 

9 = t'An-^{xi/x2). (44) 

The transformation that rotates x counter-clockwise by to a; = g{9) x onto the 
slice is found by inserting (44) into the expression for the action of SO (2) on x, 

xi = xi cos 9 — X2 sin 9 , X2 ^ xi sin 9 + X2 cos 9 (45) 
Vi = yi cos 9 — y2 sin 6* , 2/2 = 2/1 sin 9 + y2 cos 9 
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yielding the transformations in analytic form: 



X2 ~ ri = y xl + , z = z 

Vi = {x2yi- xiy2)/ri, = {xiyi+X2y2)/ri. (46) 

This transformation rotates point x into the slice point x. Alternatively, the 
transformation can be viewed as providing invariant variables on which to 
project dynamics, as we did in Hilbert basis case. Note the relation to the 
invariant polynomials (32), and observe that as the rotational degree of free- 
dom has been explicitly used, the method of moving frames requires no syzygies. 
Analytical determination of invariant variables through moving frame transfor- 
mations and particular slice conditions can be carried out systematically for 
general group representations [51], but for reasons explained in sect. 6.3, we 
shall not take this path here. 




Figure 6: State space portrait of complex Lorenz flow in reduced state space, projected on 
the slice, taking as the slice fixing point x' (a) the irreducible subspace slice condition (43), 
(b) the relative equilibrium x'q^ . 

As in sect. 6.1, note that the invariants are not well defined in the xi, a;2 — > 
limit. Using x = ri 6**^^ , y = r2 e**^^ we can write 

X2 = ri, ^ r2sin(0i - (/)2) 

^2 = r2 cos((?!)i - (/)2) , z z . (47) 

For any given y (therefore also for given 02), the limit of y for x — > docs not 
exist, as the above expression depends on the direction in the complex x-plane 
along which we approach zero. 

From a different perspective, we may redefine the slice so that xi = X2 = 
is excluded, that is by xi = 0, X2 > 0. Then we may say that group orbits of 
points in the xi = X2 = subspace fail to intersect the slice. In figure 6 (a) 
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this becomes apparent by the trajectories in reduced space being stretched as 
they come closer to xi = 2:2 = subspace, where transverse intersection would 
eventuaUy fail. 

It is instructive to rewrite the complex Lorenz equations (12) in terms of 
the invariant variables (46). This is achieved by using the chain rule (31) and 
expressing the resuh in terms of variables (46). The moving frames symmetry- 
reduced complex Lorenz equations are a 4-dimcnsional ODE system 

X2 = -(t(x2-|/2), iji = -yi+ P2X2- {e + cryi/x2)y2 

y2 = -y2 + {Pl- z)^2 + {e + (Jyi/x2)yi, Z =-6z + X2f2- (48) 

Note the singularity as 0:2 = ri — > 0. 

The projections in figure 6 help us understand the topology of the dynamics 
but also present large discontinuous jumps. Note that the invariants (46) are 
related to the invariant polynomials (32) by division by yxf + aJf. This is 
the reason we get a clearer visualization of the dynamics than with invari- 
ant polynomials: All invariants scale as the original coordinates. At the same 
time division by + x\ causes the jumps in the y components whenever the 
magnitude of x comes close to zero. 

So what docs this imply for the ultimate goal of this paper, "Continuous 
symmetry reduction and return maps for higher-dimensional flows?" In this 
example there is no problem. It is obvious by inspection of figure 6 that one can 
chose a Poincarc section away from the singular set. Repeating the construction 
of sect. 5.2 results in a return map very much like the one of figure 4 (b), with 
the same admissible orbits, so there is no need to plot further return maps for 
various choices of slices. 

6.3. Complex Lorenz equations: the general linear slice 

The irreducible subspace slice condition (43) yields the invariant variables 
(46) in explicit, analytic form. As explained in sect. 6.1, the method of moving 
frames also introduces artificial singularities in reduced state space, the location 
of which depend on the choice of slice point, and one might suspect that the 
singularities encountered by the strange attractor are due to very special choice 
of the slice fixing condition. How does a general choice of the slice fixing point 
affect the singular set? In this setting it is no longer convenient to explicitly write 
out transformations to invariant variables as we did in sect. 6; we will implement 
the moving frame map numerically, mapping computed trajectories to the slice. 
Even though the analytical computation of invariants by the method of moving 
frames can be implemented by computer algebra [20] for system dimensions of 
the order of 100, it is both computationally prohibitive and utterly unnecessary 
for symmetry reduction of very high dimensional flows of order 100, 000 required 
for fully resolved 3-D fluid simulations [57]. 

Since Qi organizes reduced space dynamics around it, but also sets the 
scale of angular velocity of symmetry induced rotations in the system, we find 
it natural to choose slice-fixing point (40) computed from group orbit of the 
relative equilibrium xq^ , given in (28) for the parameter values used here. We 
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can use (41) to compute 9 for any point x, but keeping up with the numerical 
approach of this section, we use a Newton's method. As an initial guess we use 
the angle from previous point along the trajectory, while for the first point we 
choose among the two possible solutions by demanding that the rotation brings 
the point to a minimum distance from x' . Projections of complex Lorenz flow 
to the slice defined in this manner are shown in figure 6(b). This projection 
of the strange attractor also clearly exhibits the moving frame angle 7r-jumps. 
Figure 7(a) shows the projection of singular set in 3-dimcnsions {xi,X2,yi)', the 
attractor collides with the singular set head on. 

The linear relation (42b) and figure 7 suggest that we can manipulate the 
singular set so that the attractor avoids the singularity by increasing the ratio 
yj/y'i- Choosing the slice fixing point as x' = xq^ + (5,0,0,0,0) we can 'tilt' 
the singular set so that trajectories approach it in a smoother manner, see 
figure 7 (b). 



Figure 7: State space portraits of complex Lorenz flow in reduced state space. We use a 
moving frame map to a slice fixed by point (a) x' = xq^, with the gray plane indicating the 
singular set. (b) x' = xq^ + (5, 0, 0, 0, 0). Compare with figure 6 (b). 

For more pointers on how to pick good slices and combine them with well- 
chosen Poincarc sections, the reader is referred to rcf. [20]. 

7. Differential formulation: the method of slices 

Instead of post-processing a full state space trajectory, we can proceed as 
follows: Split up the integration into a sequence of short time steps, each fol- 
lowed by a rotation of the final point such that the next segment's initial point 
is in the slice fixed by a point x' . In the infinitesimal steps limit this leads to the 
method of slices, a differential form of the method of moving frames for which 
the trajectory never leaves the reduced state space. 

Consider an A^-dimensional Lie group G acting on d-dimensional space and 
which, at least locally near x' , has ^''-dimensional orbits. For points that can be 
mapped by a moving frame to slice through x' we can write, using decomposition 



(a) 




22 



(37), the full state space trajectory as x{t) = g{T)x{T), where the (d— iV)-dim- 
ensional reduced state space trajectory x{t) is to be fixed by some condition, 
and 5(t) is then the corresponding group action on the TV-dimensional group 
manifold that rotates x into x at time r. The time derivative is i = v{gx) = gx+ 
gu, with the reduced state space velocity field given by u = dx/dr. Rewriting 
this as u = g~^v{gx) — g~^gx and using the equivariance condition (3) leads to 

u — V ~ g^^gx. 

The Lie group clement (4) and its time derivative describe the group tangent 
flow 

g-'g^g-'^e'-'^e-T. 

This is the group tangent velocity g^^gx = 9 ■ t{x) evaluated at the point x, 
i.e., with .9 = 1, see figure 5 (a). The flow in the (d—N) directions transverse to 
the group flow is now obtained by subtracting the flow along the group tangent 
direction, 

u{x) = v{x) — 9{x) ■ t{x) , u = dx/dr , (49) 

for any factorization of the flow of form x{t) = g{T)x(T). To integrate these 
equations we first have to fix a particular flow factorization by imposing con- 
ditions on x{t), and then integrate phases 9(t) on a given reduced state space 
trajectory x{t). 

Here we shall demand that the reduced state space is confined to a linear 
hyperplanc slice. Substituting (49) into the time derivative of the fixed slice 
condition (39), 

u{xft',=v{x)^t',~9a-t{xft',^Q, 

yields the equation for the group phases flow 6 for the slice fixed by x', together 
with the reduced state space M flow uilS): 

(50) 

^gAT. (51) 

Each group orbit M.x ~ {9x\g G G} is an equivalence class; the method of 
slices represents the class by its single slice intersection point x. By construction 
vFt'^ = 0, and the motion stays in the (d— A^)-dimcnsional slice. We have thus 
replaced the original dynamical system {A^, /} by a reduced system {Al, /}. 

These equations are easily integrated (provided some care is taken about how 
the trajectories cross the singular set), and given the same slice- fixing conditions, 
the integrations reproduce the plots obtained by the method of moving frames, 
such as figure 6. 

For example, consider the complex Lorenz equations slice condition of sect. 6.2: 
xi =0, X2> 0. The reduced state space equations are given by 

u{x) = v{x) - ^t{x) . (52) 



Oa{x) 

u{x) 



t{x)T ■ t' 

v{x) - 9{x) ■ t{x) , 
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Substitution of the complex Lorenz equations equations recovers (48), obtained 
by the method of moving frames. The integration of (48) recovers the strange 
attractor in figure 6 (a), obtained by simulation in the full state space, followed 
by moving frame rotation into the slice. 

In pattern recognition and 'template fitting' settings, (50) is called the re- 
construction equation. We have already encountered it in our polar coordinates 
exercise (25). Integrated together, the reduced state space trajectory (51), the 
integrated phase (50), and g(r) = exp[0(T) • T] reconstruct the full state space 
trajectory x{t) = g{T)x{T) from the reduced state space trajectory x{t), so no 
information about the flow is lost in the process of symmetry reduction. 

The denominator in (50) vanishes and the phase velocity 9{x) diverges when- 
ever the direction of group action on the reduced state space point is perpendic- 
ular to the direction of group action on the slice point x' . Therefore the method 
of slices has the same singular set as its post-processing variant, the method of 
moving frames: the intersection of the slice with the set of points with group 
tangent perpendicular to t' . 

A historical note. The basic idea of the method of slices is intuitive and 
frequently reinvented, often under a different name; for example, it is stated 
without attribution as problem 1. of Sect. 6.2 of Arnol'd Ordinary Differen- 
tial Equations [58]. The factorization (37) is stated on p. 31 of Anosov and 
Arnol'd [59], who note, without further elaboration, that in the vicinity of a 
point that is not fixed by the group one can reduce the order of a system of 
differential equations by the dimension of the group. Fiedler, in the influential 
1995 talk at the Newton Institute, and Fiedler, Sandstede, Wulff, Turaev and 
Scheel [60-63] treat Euclidean symmetry bifurcations in the context of spiral 
wave formation. The central idea is to utilize the semidirect product structure 
of the Euclidean group E{2) to transform the flow into a 'skew product' form, 
with a part orthogonal to the group orbit, and the other part within it, as in 
(51). They refer to a linear slice near a relative equilibrium as a Palais 
slice, with Palais coordinates. As the choice of the slice is arbitrary, these co- 
ordinates are not unique. According to these authors, the skew product flow 
was first written down by Mielke [64], in the context of buckling in elasticity 
theory. However, this decomposition is no doubt much older. For example, it 
was used by Krupa [28, 45] in his local slice study of bifurcations of relative 
equilibria. Biktashev, Holden, and Nikolaev [65] cite Anosov and Arnol'd [59] 
for the 'well-known' factorization (37) and write down the slice flow equations 
(51). Haller and Mezic [66] reduce symmetries of three-dimensional volume 
preserving flows and reinvent method of moving frames, under the name 'or- 
bit projection map.' There is extensive literature on reduction of symplectic 
manifolds with symmetry; the Marsden and Weinstein 1974 article [67] is an 
important early reference. Then there are studies of the reduced phase spaces 
for vortices moving on a sphere such as rcf. [68], and many, many others. 

Neither Fiedler et al. [60] nor Biktashev et al. [65] implemented their meth- 
ods numerically. That was done by Rowley and Marsden for the Kuramoto- 
Sivashinsky [53] and the Burgers [69] equations, and Beyn and Thiimmler [70, 
71] for a number of reaction-diffusion systems, described by parabolic partial 
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differential equations on unbounded domains. We recommend the Barkley pa- 
per [72] for a clear explanation of how the Euclidean symmetry leads to spirals, 
and the Beyn and Thiimmler paper [70] for inspirational concrete examples of 
how freezing/slicing simplifies the dynamics of rotational, traveling and spiraling 
relative equilibria. 

Beyn and Thiimmler write the solution as a composition of the action of a 
time dependent group element g{T) with a 'frozen,' in-slice solution u{t) (37). 
In their nomenclature, making a relative equilibrium stationary by going to a 
co-moving frame is 'freezing' the traveling wave, and the imposition of the phase 
condition (i.e., slice condition (38)) is the 'freezing ansatz.' They find it more 
convenient to make use of the equivariance by extending the state space rather 
than reducing it, by adding an additional parameter and a phase condition. 
The freezing ansatz [70] is identical to the Rowley and Marsden [69] and our 
slicing, except that freezing is formulated as an additional constraint, just as 
when we compute periodic orbits of ODEs we add Poincare section as an addi- 
tional constraint, i.e., increase the dimensionality of the problem by 1 for every 
continuous symmetry. 

Our derivation of method of slices follows most closely Rowley and Mars- 
den [69] who in the pattern recognition setting refer to the slice point as a 
template^ and call (50) the reconstruction equation [73]. They also describe the 
'method of connections' (called 'orthogonality of time and group orbit at succes- 
sive times' in ref. [70]), for which the reconstruction equation (50) denominator 
is i(a;)"^ •t(x) and thus nonvanishing as long as the action of the group is regular. 
This avoids the spurious slice singularities, but it is not clear what the method 
of connections buys us otherwise. It does not reduce the dimensionality of the 
state space, and it accrues geometric phases which prevent relative periodic 
orbits from closing into periodic orbits. Geometric phase in laser equations, 
including complex Lorenz equations, has been studied in ref. [74-78]. 

One would think that with all the literature on desymetrization the case is 
shut and closed, but not so. Applied mathematicians are inordinately fond of 
bifurcations, and almost all of the previous work focuses on equilibria, relative 
equilibria (traveling waves), and their bifurcations, and for these problems a 
single slice works well. Only when one tries to describe the totality of chaotic 
orbits does the non-global nature of slices become a serious nuisance. 

8. Discussion and conclusions 

We have presented two approaches to continuous symmetry reduction of 
higher-dimensional flows and illustrated them with reductions of the complex 
Lorenz system, a 5-dimensional dissipative flow with rotational symmetry. In 
either approach numerical computations can be performed in the original, full 
state-space representation, and then the solutions can be projected onto the 
symmetry-reduced state space. 

In the Hubert polynomial basis approach, one transforms the equivariant 
state space coordinates into invariant coordinates by a nonlinear coordinate 
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transformation 

{xi,X2, ■■ ■ ,Xd} ^ {ui,U2, ■ ■ ■ ,Um} , 

and studies the invariant image of dynamics rewritten in terms of invariant co- 
ordinates. These invariant polynomial bases can be algorithmically determined 
for both Hamiltonian and dissipative systems. Our goal is to reduce symmetries 
of fully-resolved simulations of PDEs, with state space dimensions of the order 
of a few tens to few hundreds (for Kuramoto-Sivashinsky flow), and well into 
the tens or hundreds of thousands (for pipe and plane Couette flows). Unfor- 
tunately, the computational cost of polynomial basis algorithms is at present 
prohibitive for state space dimensions larger than ten, so the invariant polyno- 
mial basis approach is not a feasible option. We have discussed it here solely 
for illustrative purposes. 

In the method of moving frames (or its continuous time, differential version, 
the method of slices), one fixes a local slice ix — x'^t' = 0, a hyperplane normal 
to the group tangent t' that cuts across group orbits in the neighborhood of the 
slice- fixing point x' . The state space is sliced locally in such a way that each 
group orbit of symmetry-equivalent points is represented by a single point, with 
the symmetry-reduced dynamics in the reduced state space M/G given by (51): 

u = v-e-t, e^{v-t')/{t-t'). 

The method of moving frames turns out to be an efficient method for reducing 
the fiow to a symmetry-invariant reduced state space, suited to reduction of 
even very high-dimensional dissipative flows to local return maps: one runs the 
dynamics in the full state space and post-processes the trajectory by the method 
of moving frames. Importantly, from a numerical point of view, there is no need 
to actually recast the dynamics in the new coordinates or write new code. Either 
approach can be used as a visualization tool, with all computations carried out 
in the original coordinates, and then, when done, projecting the solutions onto 
the symmetry reduced state space by post-processing the data. In contrast to 
co-moving frames local to each traveling solution, restricting the dynamics to a 
slice renders all relative equilibria stationary in the same set of coordinates. 

An inconvenience inherent in the linear slices formulation is that they are 
local, and the reduced flow encounters singularities in subsets of the reduced 
state space, with the reduced trajectory exhibiting large, slice-induced jumps. 
This singular set is introduced by and depends on the slice-flxing condition. We 
have shown, in the 5-dimensional complex Lorenz equations example, that the 
location of the singular set can be manipulated by judicious choice of the slice 
fixing point, and geometrical information about the dynamics can be extracted 
by constructing a return map through a Poincare section that does not inter- 
sect the singular set. The trick is to construct a good set of symmetry invariant 
Poincare sections, and that is a dark art for systems of dimension higher than 
three, with or without a symmetry. In higher-dimensional fiows, with more in- 
volved symmetry group actions and larger sets of stationary solutions, where 
a single slice and Poincare section will not suffice, we can still expect to cover 
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the reduced state space with multiple slices, obtaining a set of discrete maps in- 
volving multiple Poincare sections. As to high-dimensional applications, it was 
shown in ref. [20] that the coexistence of four equilibria, two relative equilib- 
ria and a nested fixed-point subspace structure in an effectively 8-dimensional 
Kuramoto-Sivashinsky system [19] complicates matters considerably. This ap- 
plication of symmetry reduction to a spatially extended, PDE system is the 
subject of a forthcoming publication [79]. 
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